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1. INTRODUCTION 



In superstring as in grand unification, the high energy extrapolation of the standard 
model renormalized gauge coupling constants is described by a one loop scale evolution of 
familiar form: 

|g = (^ +2()alog £_ + A„(M i ,M i ). (!) 

(The index a = 3,2,1 labels the SU(3) x SU(2) x U(l) group factors G a and b a are 
the beta function slope parameters associated with the low energy modes, f3 a (g) = 
— 6 a (7^/(4-7r) 2 + • • •) The superstring case is, however, distinguished by three important 
features [1]: 

(i) Tree level relations [2] involving the gauge and gravitational interactions, 

2 2 _ 4K 2 __ 327T 

9x - k a9a - - ^mJ- (2) 

In addition to the string theory expansion parameter gx (or 4-d dilaton VEV < S >= 
l/g x ) which is specified by the ratio of the string mass scale Ms = -j= to the phenomeno- 

logical Planck mass M P = ^ = 1.22 x 10 19 Gey, as exhibited in eq.(2), three extra free 
(positive rational numbers) parameters k a are introduced in eq.(2), corresponding to the 
levels of the affine Lie algebras for the gauge group factors G a in the underlying string 
theory. 

(ii) An improved unification scale Mx defined in eq.(l) as the matching scale between 
the field and string theories renormalized coupling constants at which these obey most 
closely the tree level relations, eq.(2). For the field theory coupling constants in the DR 
regularization scheme, 

e (l-7)/2 e (i-7)/2 
M x = ——=g x Mp = ——Ms ~ Sx5.27 x \Q 17 GeV. (3) 
47TV27 yZTiyZl 

The field theory (ft) convention in use here is related to the string theory (st) one as, 
g[ l = \Plg S a \ corresponding to the normalization of the Lie algebra generators, Trji(Q 2 ) = 
^c(R), where c(R) = l(R) is the Dynkin index of representation R. 

(iii) Threshold corrections accounting for the contributions of the infinite set of mas- 
sive string states at the string (Ms) and compactification (Mc) scales, integrated out 
by matching the field and string theories scattering amplitudes. These corrections are 
represented in eq.(l) by the functions A a (Mi, Mi) depending upon the structure of the 
string mass spectrum and the other characteristic parameters of the compactified space 
manifold, such as the vacuum expectation values (VEVs) of the compactification moduli 
fields, Mi = Tj, Ui [3]. Specifically, Mx is defined as the choice of scale which minimizes 
the threshold corrections contributions. Of course, the perturbative character of formula 
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(1) implies that the size of A a should be comparable to that of two- loop effects, so that 
A a = 0(l). 

For a quantitative test of superstring unification based on eq.(l) and for a proper 
identification of the fundamental parameters, Mx and gx, it is essential to understand 
the structure and size of threshold corrections. Thus, an additive decomposition such as 
A a = k a Y — b a A, may be exploited to introduce effective unification scale and coupling 
constant, 

M x ^M' x = M x e N \ 9x^9'x = ( 4 ) 

so defined as to incorporate the contributions from the above two components, Y and A. 

The toroidal compactification orbifold models prove very helpful in obtaining an in- 
formation on A a . The contributions from compactification modes admit here a natural 
additive decomposition into a moduli dependent component arising from the chiral mass 
F-terms and a moduli independent component arising from the vector mass D-terms [4]. 
As is well known, the moduli dependent contributions play an essential role in the cancel- 
lation of sigma model anomalies affecting the target space duality symmetry [5]. These can 
be represented by general formulas involving the automorphic functions of the compacti- 
cation manifold accompanied by model dependent coefficients. On the other hand, the 
moduli independent contributions carry only an implicit dependence on the compactifica- 
tion manifold, such as the orbifold gauge embedding or the discrete Wilson lines. In spite 
of several attempts in the literature to estimate the size of both components of threshold 
corrections [1,6-9], one is still lacking a clear physical understanding of their magnitude. 
Our main goal in this paper is to present results for the moduli independent threshold 
corrections through an extensive numerical study based on a sample of orbifold models. 
A recent work by Dienes and Faraggi [34], which appeared while this paper was being 
completed, pursues a similar goal to ours based on the fermionic models. 

The main physical motivation for this paper is, however, the wide gap that separates 
the improved string unification scale Mx = O.2I6M5 ~ 5 x lQ 17 GeV, assuming gx = 
O(l), from the observed grand unification scale, Mqut - 2 x l^GeV, as determined 
by extrapolating the gauge coupling constants up from their experimentally determined 
values at the Z— boson mass [10]. The implications of this order of magnitude discrepancy 
in scales have been emphasized on several occasions [11]. The conflict for superstring 
unification can be resolved in two different ways: One can postulate [12] large string 
threshold corrections such that after becoming equal and joining together at the observed 
scale Mqut the gauge coupling constants follow diverging flows up to Mx- A matching 
of the one loop extrapolated values of g a (Mx) with their predicted values, as obtained 
by adjusting the moduli dependent threshold corrections, can be successfully achieved in 
terms of wide classes of solutions for the modular weights of massless modes consistent 
with the anomaly cancellation constraints [5,12]. Alternatively, one can postulate [13] an 
affine level parameter for the weak hyper charge group U(l)y somewhat lower than the 
standard grand unification group value, k\ = |. With such an enhanced starting value for 
{kia\ (mz)) -1 one achieves a delayed joining of the gauge coupling constants flows which 
can easily raise up the unification scale by one order of magnitude. While either of these 
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possibilities is well motivated by itself and appears sufficient to rescue a superstring grand 
desert scenario, there remains certain unsatisfactory points. Thus, the VEVs of moduli 
fields requested in the first possibility, < T >= 10 — 30, appear to be somewhat too large, 
whereas no known realistic orbifold example [14,34] exists for which the hypercharge group 
level parameter comes as low as the value k\ ~ 1.4 favored in the second possibility. 

A generic feature of standard-like orbifold models is the occurrence of a rich spectrum 
of charged massless modes appearing on side of the requested (quarks and leptons) chiral 
families in vector representations of the color and weak groups. In fact, the matter rep- 
resentations of the observable sectors group factors are generally sizeable enough so that 
the corresponding beta function parameters (3 a arise with either small negative values 
or large positive values. This suggests that a first stage of slow or non asymptotically 
free scale evolution may well take place from Mx down to some scale where the extra 
modes pair up by acquiring mass and decouple. As is well-known [15], in order for the 
4-d low energy effective theory to be weakly coupled, so as not to invalidate the use of 
eq.(l) (gx ~ gaM^ < 1, ga = 10-dimensional gauge coupling constant), and in order to 
avoid dealing with a strongly coupled 10-d theory (gdMg < 1) one must require that the 
compactification and unification scales retain a magnitude comparable to the string scale, 
Mx — Mq — Ms- (The second restriction can be relaxed by allowing, for instance, for an 
anisotropic compactification manifold (large radius in one out of the six compactified di- 
mensions) in which a weakly coupled effective theory, gx < 0(1), could remain compatible 
with a strongly coupled string theory (large go) [16].) Assuming the above near equality 
of scales, a natural identification for the decoupling scale of the extra matter is the mass 
scale, denoted Ma, which is induced by a non- vanishing Fayet-Iliopoulos D-term contribu- 
tion to some apparently anomalous £7(1) group factor occurring on compactification [17]. 
This suggestion is not new, of course, and appears in several places in the specialized liter- 
ature. The idea is to cancel the non vanishing one loop string contributions to the D-term 
scalar potential of an apparently anomalous U(l) factor by judiciously lifting the VEVs of 
certain scalar fields while restoring a stable supersymmetric vacuum. We shall carry out 
an analysis of the one loop gauge coupling constants unification which combines together 
the above ideas of adjustable moduli VEVs and k\ level parameters together with that of 
an adjustable intermediate scale Ma, while describing the scale evolution in the interval 
from Mx to Ma on the basis of orbifold models predictions. 

The paper contains 5 sections. In Section 2, we discuss in wide outline the basic 
formalism involved in the one loop string renormalization of the gauge coupling constants 
as applied to orbifold models. None of the results discussed in this section is new, our 
main intent being to provide a concrete, encapsulated presentation of the relevant formal- 
ism. In Section 3, we present numerical results for the moduli independent threshold 
corrections for a sample of representative orbifold models. In Section 4, we examine the 
viability of superstring unification in an extended picture including threshold corrections 
and an intermediate scale associated with an anomalous U(l) symmetry. In Section 5, we 
summarize the main conclusions. 
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2. ONE LOOP STRING RENORMALIZATION 



2.1 Threshold corrections to gauge coupling constants 



We consider the class of low energy supersymmetric theories descending from 4-d 
heterotic string theories with a nonsemi-simple gauge group G a - The genus zero (unity) 
world sheet (with Wick-rotated Euclidean metric) of the conformal field theory is a sphere 
(torus) parametrized by planar coordinates: with corresponding 

cylindrical coordinates given for the sphere by: £ = a — it, £ = a + it, a G [0, 1], t G 
[— oo, oo] and for the torus by: £ = a + rt, C = a + ft, a, t G [0, 1], r = T\ + ir 2 . The 
right-moving RNS (Ramond-Neveu-Schwarz) superstring is built with 20 spacetime and 
spin fields X^(z), ip^(z), [ji = 0, •■■,9], associated with D = 4 external dimensions of 
the flat spacetime [/j, = 0, • • • , 3] and d— D = 10 — D = 6 internal dimensions [/j, = 4, • • • 9] 
of the compactification space manifold, represented in a complex basis as, X R , X l R , ip l = 
e^^ip 1 = e -1 ^, [i = 1,2,3] where the complex scalar fields <j>i(z) are coordinates of the 
SO (6) group Cartan torus. This is tensored by a left-moving bosonic string built with 26 
fields X^(z), [fi = 0, 25], comprising D external space coordinates and 26 — D internal 
space coordinates which are distributed into 6 compactified space coordinates X^^X"^ and 
16 gauge coordinates of the E 8 x E' 8 Cartan torus F^F ^i" = 1,...,8] generating the 
currents J a (z) of the affine Lie algebras G a of levels k a . At certain places, we refer to these 
coordinates globally as F 1 , [I = 1, • • • , 16] and also by using their fermionic representation 
in terms of complex 2-d Weyl spinors, (A a ,A Q ) = e ±lF , [I = 1, - - • , 16; a = l,---8]. 
Of course, the above covariantly quantized string theory must be supplemented with the 
conformal ghost fields c z (z, z), b zz (z, z) and the superconformal ghost fields, j(z),P z (z) 
[18]. 

The one loop string threshold corrections are described by a general formula obtained 
by Kaplunovsky [1], 

A a = k a Y + A a , A a = - JjLL^k a B a {q,q) -b a ^j, (5) 

where one has decomposed the total contribution, denoted A a , into a universal contribu- 
tion, k a Yo, independent of the gauge group factor (except for the coefficient k a ), arising 
from gravitational interactions and oscillator excitations modes, and a contribution solely 
due to the massive compactification modes, denoted A a . The latter component is expressed 
as a deformed partition function integrated over the inequivalent complex structures of the 
genus 1 world sheet, with an integrand 



d (T) \ 
2 ^ I ?7(r) 2 ?7(T) 2 dq\ t]{t) J 

even ct,p 



(6) 



x 2Trace(^(-l) 2 ^ F Q 2 a q Lo -^q L "-^ : 
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where the first factor represents the partition function of the external theory inserted with 
the operator ( — + x 2 ), where \ denotes the 4-d helicity or chirality vertex operator and 
we have introduced the familiar Dedekind function, i](r) = J| n (l — q n ), and the Jacobi 
theta-functions (cf. eq.(10) below). 

The second factor in eq.(6) (with F = fermion number operator, Lq,Lq = conformal 
dimensions operators) corresponds to the internal theory partition function inserted with 
the square Q\ of any one of the gauge group generators for subgroup G a . The integral over 
the world sheet torus complex modular parameter, r = ti+zt2, with q = e 27riT , q = e _27r * T , 
extends over the modular group SL(2,Z) fundamental domain, F = [\ti\ < |,|r| > 
1]. Infrared convergence of the integral, eq.(5), is ensured by the subtraction of b a = 
lim^^oo k a B a , where b a = \ E a h c s(^a) - 2c F ( J R Q ) + llc v (R a )] (S = complex scalar, 
F = Weyl or Majorana fermion, V = vector) represent the summed contributions to the 
beta function slope parameters from the massless string modes a. 

The summation in eq.(6) over the subset of even spin structures of the right- moving 
sector, (a,j3) = [(0, 0), (0, §), (±, 0)] where a,j3 = = NS(A) (Neveu-Schwarz, Antiperi- 
odic) or | = R(P) (Ramond, Periodic) is performed by insertion of the familiar GSO 
(Gliozzi-Scherk-Olive) projection phase factors leading to the supersymmetric string [18]. 

2.2 Specialization to orbifolds 

To express the second internal space factor in eq.(6) for orbifolds, we recall first that 
the projection (modding) with respect to the orbifold point symmetry is achieved by 
summing over the (space and time) twisted subsectors (g, h) by using [19,20], 

Trace(- ■ •) = ]T ]T h)Trace g (h ■ ■ •), 

g h;[g,h]=0 

where |G| is the orbifold point group order and x(g, h) are degeneracy factors. For toroidal 
compactication, all fields are free so that the torus partition function is obtained by as- 
sociating to a complex coordinate field X(a,t) of given chirality, a factor 1/ (tt\/2t2) (flat 
case) or (1 — e 27riv ) /i](t) (untwisted case with time twist X(a,t + 1) = e 27rtv X(a, t) ) or 

v{ T )/$[l + y] (space twisted case X(a+ l,t) = e 27Zlv X(a,t)) and to a fermionic Majorana- 

Weyl field, obeying the twisted boundary conditions: 

V>(<7 + l,t) = -e 27 " 9 V(a, t), *P(a, t + 1) = -e" 27r ^V((7, t), 

a factor ["#[^/]/?7(t)]5. The zero modes are associated a factor qPL/ 2 qPn/ 2 summed over 
the winding modes spanning the compact ificat ion manifold lattice Aq with basis vectors 
e\ and over the Kaluza-Klein momentum modes spanning its dual lattice A* with basis 
vectors e* a (cf. eq.(ll) below)). 

We recall next that a torus _R 6 /A 6 , defined by X 1 = X 1 + 2nn a e l a: having a point 
symmetry group, P = Zjv, of automorphisms of the lattice A 6 , defines an abelian orbifold 
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endowed with a space symmetry group, G = P x A 6 . The space group action on the string 
theory fields is described in terms of rotations 9 k and translations u k j together with their 
associated gauge group shift embedding elements described by translations V 1 and Wilson 
lines translations a^, [I = l,---,16;a = 1,2,3]. The space group G = {g p } = {f3 p ,w p } 
composition laws read: g x g 2 = P1W2 + ^i), 9 P 1 = (Pp 1 , -fi^Wp). 

The string Hilbert space of states consists of the untwisted sector (k = 0) and the 
twisted (k = 1, • • • , N—l) sectors. The twisted sectors g k are distinguished by the boundary 
conditions: (X(a + 1, t), i/)(a + 1, t)) = g k (X(a, t), -(-l) 2 ^(a, t)). They are organized 
into conjugacy classes of the space group with representative elements, g k = [9 k , u k j] and 
their associated classes, {g k ~ g'g k g _1 = (9 k ,u k ),g' = g p G Zn}, where the set of shift 
vectors u k = (6 p Ukj + (1 — 9 k )u), [u G A 6 ,p = 0, • • • N — 1], span lattice cosets (labelled 
by the index /) with representative elements, Ukj- The compactified space coordinates, 

X 1 = X L +X l R = x 1 +i7ttp l + 27raw l H (units 2a' = 1), admit the (zero and oscillators) 

modes expansion, 

mi 

In twisted sectors, the string center of mass coordinates x l are not arbitrary real parameters 
but rather must satisfy: g k x = x + Ukj + u, [ukj,u G A 6 ]. Therefore, each of the g h 
twisted sectors splits into subsets which can be classified in terms of the corresponding set 
of fixed points of the space group, f^ 1 , defined as: 9 k f^ = /( fc ) + u k j where u l k * = 

m t f e ai [ m<1 = integers] are translation vectors of the 6-d toroidal lattice Aq determined 
by the condition that they return the rotated fixed point 9 k f back to its original position, 
so that / = (1 — k )~ 1 Ukj + u. Specifically, the k-twisted sector fixed points are 
distinguished by a label a running over the number of fixed points. The lattice vectors 
Ukj identify with the lattice coset representatives u k j introduced above only for prime 
orbifolds. For simply twisted sectors, k = l or k = N — 1 = —l(modN), the fixed points 
and conjugacy classes u k j are in one to one correspondence, so that faithfully 
label these classes and u k j = u k j. This property holds true for all the twisted sectors 
in the prime orbifolds, Z 3 j. For the multiply twisted sectors, the full set of fixed points 

decomposes into disjoint subsets {f^ \ f} k \ ' ' •}> where the fixed points within each 

subset (labelled by A) are related as, 9 PA f^ = fj^ 7^ for pa < k, and hence are in 
one to one correspondence with the same conjugacy classes, u k j. The cases involving non 

trivial subsets [f^], comprising more than one fixed point, arise only for the non prime 
(N = 1) orbifolds .£4,6,8,12 and for the direct product orbifolds Zn x Zm- 

The orbifold space group elements can now be expressed as , 

g k = {9 k , u kJ = m a kJ e a - W 1 = kV 1 + m a kJ a? a }, 



9 k = diag(9*) = diag(e 2lTlkv *) 



5> = ° 



(7) 



The orbifold group action on fields (eq.(8)) and state vectors (eq.(8')) reads in obvious 
notations: 

g k X* LiR = 9{X i LiR + 2irml f el g k F I = F 1 + 2*^* + m^oj), g% = 0,ty<, (8) 

/[(«i ni )^[(«'-m,) 9i ](-)bi? 5 ^ = « l + >i? Ipl.P 7 = + ^ >L 

The above used correspondence between Wilson lines translation vectors and the non con- 
tractible loops, refers to abelian orbifolds. Non abelian orbifolds with shift gauge 

embeddings can be constructed by extending the definition of Wilson lines to class depen- 
dent shift vectors, kV 1 — > V/j derived from a gauge embedding matrix of general form 
[21]. 

The internal space oscillator operators, (a l n . , a J mj )( L y where i, j are complex conjugate 

bases indices, (given by the familiar linear combinations of real basis indices, \i = (1 + 
i2)/v / 2, (1 — i2)/y/2, ■ ■ •), enter with the moddings, rii G Z =F Oi, rrij E Z ± 9j, where Z 
designates the set of integers. The translation vectors a 1 = (rii + |), [rii G £V m G 
2Z+l(odd integers)] are elements of the 50(6) group weight lattice T 6 and W 1 = n 1 , (n 7 + 
^), [n 7 G Z 1 Y^i=\ nI £ (even integers) ] are elements of the Eg x -Eg group weight 
lattice, Tg+g. The translation vectors v l and V r ,a^ with respect to these lattices must 
obey: Nv l G r6,iV"l/ 7 G Fs+s, Nm a a^ G Tg+g as well as the level matching (modular 
mvanance under T N ) conditions N^kV 1 + m^afj 2 - (kv 1 ) 2 } G 2Z. 

With the above rules in hand, we can now quote the following more explicit formula 
derived from eq.(6): 



B a(g,?) = -2|^5^x(TO»rc)e(TO,n)- 



* n 

i=l,3 



1 1 

X - 



euena,/? 



("I) 



2a+2/3 



-2g- 



r 







/3+nt>i 



?/(r) 



n 

i=l,3 



7/(r) 



(f)0 



(r) 



x 



4t/ 16 (t) 

n ^ 



r](m,n;a,f3;a',/3') JJ Q 72 ^ 



J=l 



a' + mVjf 



j=i 



P£/2gP?t/2 



/3 + nV> 



(r) 



■i] 2 (t)i] 2 (t) dq \ r)(r) 
t/(t) 



(9) 

where the second and third factors, recognizable by the brackets, are contributed by the 
internal space coordinates and spinors, the fourth factor by the gauge coordinates and the 
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last (fifth) factor by the compactified space zero modes. The numerical factors appearing 
in denominators account for the averaging over the time-like spin structures. 

2.3 Classification of threshold corrections 

The generalized GSO orbifold projection, which selects the singlet states with respect 
to the orbifold space symmetry group, is represented by the sum over the various twisted 
orbifold subsectors, (g, h) = (m, n), performed jointly with the sum over the spin structures 
(a, /?), (a/, /?') for the fermionized fields associated with the gauge degrees of freedom. 

The summations over twisted subsectors (m, n), (a, (3), (a/, (3') are weighted by phase 
factors e(m, n) and ?y(m, n; a, (3; a', (3') which are determined by the requirement that T 2 B a 
be invariant under the modular SL(2,Z) group, generated by S : r — > — - and T: r — > 
r + 1. The set of twisted (g, h) subsectors are mixed together under the action of the 
modular group according to the transformation law [19,22]: r — > (ar + b)/(cT + d); (g, h) — > 
(h c g d , h a g b ), [a, 6, c, d G Z, ad - 6c = 1]. (For Z/v orbifolds, 5 : (m, n) -> (N - n,m), T : 
(m, n) — > (m, m + n).) The entire set of twisted subsectors can be organized into disjoint 
subsets (orbits) of subsectors which close under the modular group action. The inter-orbit 
phase factors rj(m, n, ...) are fixed uniquely by the requirement of modular invariance. The 
intra-orbits (discrete torsion) phase factors e(m, n) are independently fixed by constraints 
derived from higher string loops modular invariance or from unitarity [23] . The additional 
freedom that might be present when the factors e(m, n) are non-trivial phases serves then 
to label distinct string theories constructed from the same orbifold. Orbifolds with no (g, h) 
fixed 2-d torus (i.e., not simultaneously fixed by both space g and time h twists) possess 
one modular orbit only. Orbifolds having one simultaneous (g, h) fixed 2-d torus possess 
several modular orbits which are in correspondence with the distinct N = 2 suborbifolds 
of the initial orbifold. 

The multiplicity factors x(m, n) = x(g, h) count, for twisted subsectors, the number of 
distinct degenerate subsectors associated with fixed points of the orbifold point group which 
are simultaneously invariant under both g and h [24]. (Useful information on these factors 
is provided in refs. [25,26]). For untwisted sectors (m = 0), there occurs corresponding non 
trivial factors x(±> h) from the projection on oscillator states symmetric with respect to 
the orbifold point group. These can be explicitly calculated from the formula: x(l, 6 n ) = 
Y[i | — 2z sin(7rm>i)| 2 = |det'(l — 9 n )\, where the product and determinant are understood 
to extend over the rotated 2-d tori planes. 

In the presence of Wilson lines, an additional summation must be included over the 
independent Wilson lines a a satisfying the property 9 k a a ^ a a and over the independent 
noncontractible loop parameters labeled by m a . The overall sum over twisted subsectors 
in eq.(9) is then replaced as, J2 m n = J2 a Ylmnm a • F° r the abelian direct products 

' aii m,f 

orbifolds, Zn x Zm, [M = pN,p e Z] straightforward extensions of the above rules ap- 
ply in which one deals with pairs of generators, (#1,02)5 shift vectors, (y±, V2), (Vi, V2), 
twisted subsectors, {g\gi\ ^1^2) = ("7.1TO2; ^177,2), setting the discrete torsion phase factor 
as [14,23], e(miTO 2 , nm 2 ) = e^Hm^-m^/N ? [k = 0, ■ ■ ■ , N - 1]. 
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The inter-orbit phases r}(m, n; a, (3; a', (3') depend, of course, on the conventions adopted 
for the fermionic determinants. The following carefully chosen phase conventions for theta- 
functions [23,27], 



detd, 



a/3 



(v = 0|r) = e 



a + 9 

P + <P 



(u = 0\t), 



(10) 



nez 



which describes the determinant of a free complex Weyl field obeying the boundary con- 
ditions specified a few paragraphs above, is found to reduce the modular invariance con- 
straints on the coefficients to the remarkably simple solution of unit phases, rj(m, n; ...) = 1. 
To prove this statement in the orbifold case, one can follow the same steps as in [27] 
involving the use of the identities relating the fermionic and bosonic representations of 
theta-functions and of the Poisson formula transforming the summation over the compact- 
ification lattice to that over its dual. Combining in this way the fourth (gauge sector) 
and fifth (zero modes) factors in eq.(9) yields an equivalent representation for the product 
of these factors in terms of a manifestly modular invariant sum over an even, self-dual 
(shifted) (22,6)-dimensional Lorentzian lattice, 



weA 6 ,peA*,wer 8+8 



Pl,r = [pl m , Pi = W T + kVr,p Rli ], 



p L > R = ±G^w l 



+ - 



2B^w l 



1 



(11) 



where w. 



Pr) = m a e^,p^ = p Lil +p Rtl = n a el a , [m a , 



n. 



winding and 



momentum modes integers], G^ X G\ U = 5^ and the basis vectors norms Ylfi( e a) 2 identify 
with the compactification radii R a . The background metric and antisymmetric tensor 
fields, (G M j, , B = (G a 6, S a 6)e* a e* 6 , and the Wilson line vector field, A 1 ^ = a^e* a , repre- 
sent the generalized coupling constants of the world sheet sigma model of the heterotic 
string whose action (specialized to the superconformal gauge) is reproduced below, for 
definiteness, 

J^7 / / d(Tdt y/hh a ^d a X^X^ + i^ R p a V^G^(X) 



+6 a P(d a X> J <dpX»B^(X) + d^dpFxAUX) ) - a'VhR^D(X) 



(12) 



where V ' a ip u = d a ip v + Q^d a X x ipL l ', [O = generalized spin connection with respect to to 
the metric and torsion tensors] and D(X) = —^lnS(X) denotes the dilaton field. The 



10 



a— model background fields in orbifolds, as in toroidal manifolds, are X— independent 
constants, due to the vanishing curvature tensor. 

The charge generators Q a in eq.(6) identify with the zero modes components of the Lie 
algebra G a gauge current vertex operators, Q a = J° = J ^jJ a (z). The allowed currents 

are chosen among the linear combinations of the vertex operators, {idF 1 (z), e lPlF 
invariant under the orbifold group. Any choice of component Q^W = !>•••> dim(G a )], is 
admissible since all the components squared Q% 2 contribute equally to the trace over string 
states. It is easiest to work with the Cartan subalgebra generators because of the simpler 
structure of their representation as linear combinations of the momentum operators, Q a = 
Qai I t^ idFl with coefficients Q aI such that Q ai = Q a iE{, (£/, E}\ [i = 1, • • • , 16] 
are the moving orthogonal frames basis and its dual for the Tg+s torus) represent the 
directions (flat components) in the E 8 x E' 8 weights lattice invariant with respect to the 
orbifold group subject to the constraints, QaiV 1 \Q a io-i £ Z. The weight lattice vector 
components representing the eigenvalues of the Cartan subalgebra operators, Q", [a = 
1, • • • , rank(G)], for the momentum eigenstates, \P I = W 1 + kV 1 >, are given by the 
scalar products: {Q^-P = Q^jP 1 }. These relations can be used to explicitly determine the 
Qai, their absolute normalization being fixed by reference to the normalization condition, 
Tr(QaQb) = \c{R)5 a b, for the associated matrices. 

For non-abelian subgroup factors, the gauge group shift embedding case, to which 
we have limited our considerations here, always leads to k a = 1. For abelian subgroups, 
the parameters k a , which are still called levels for convenience of language, depend on 
the normalization of the corresponding charge operators Q a and specified by [14]: k a = 

The insertion of the charge squared operators is accounted for, in the notations intro- 
duced in eq.(9), by replacing the theta-function factors by modified ones using the following 
rule: 

8 8 

W^iQi 2 ^ E QiQ J M II ^ + E(Qa) 2 tf/'n^ (13) 

/ I^J=1 K^I,J 1=1 K^I 

where the primed and double-primed theta-functions are defined in terms of the sum 
representation given in eq.(10) by inserting linear and quadratic powers of the lattice 
momenta according to the prescriptions: 

$ = Y^q p2 f\ 9' 1 = ^?V 2/2 , 9" I = 2q^ 1 = Y,P I2 q p2/ \ (14) 
p p Q P 

using self-evident shorthand notations. Note that the precise definition of the 4-d chirality 
operator, introduced after equation (6), reads in these notations, 



9 d , •& 1 It? no 7 



n=l 



The rules in eqs.(13) and (14) follow directly from a consideration of the bosonic repre- 
sentation of the partition function, as described above in connection with eq.(ll). We 
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caution, however, that these rules are not sufficient by themselves in dealing with cases 
involving massless charged oscillator states. For these fortunately rare cases, one needs 
to insert proper correction factors in order to ensure a correct normalization of the beta 
function parameters. 

Turning now to the threshold corrections as calculated from eq.(9) we note that the 
A a have a natural additive decomposition in terms of moduli dependent and independent 
contributions which we associate to the first and second terms in the formula: 

A a (M,M) = 8 a + A^(M,M). 

This separation arises when one classifies contributions according to the number N = 4, 2, 1 
of space-time supersymmetries which are realized in terms of disjoint subspaces of the 
Hilbert space of states [3]. There exists a one to one correspondence between the super- 
symmetry irreducible representation spaces and the spaces of states of suborbifolds which 
are constructed from subgroups of the full point symmetry group, themselves identified 
with the modular orbits. The N = 4, 2, 1 supersymmetries are then associated with the 
suborbifolds leaving fixed 3, 1 or 2-d tori, respectively. The N = 4 supersymmetric 
subsector arises from the purely toroidal, trivial orbit, (g, h) = (1, 1), which is clearly ab- 
sent in orbifolds, due to the projection. The moduli dependent terms originate from N=2 
suborbifolds (one fixed 2-d torus) subsectors and the moduli independent ones from the 
N=l suborbifolds (no fixed 2-d torus) subsectors. The N = 1, 2 orbits generally contribute 
to both b a or A a while the N = 4 toroidal subsector (g, h) = (1, 1) (three fixed 2-d tori) 
contributes to neither. 

The moduli dependent N = 2 contributions arise necessarily from subsectors having 
non-vanishing momenta, Pl,r- Indeed, a non-trivial zero modes factor different from unity 
occurs only for twisted subsectors (m, n) with a simultaneous fixed 2-d torus. For this case, 
the factors in the partition function in eq.(9) multiplying the zero modes factor combine 
into the product of an holomorphic function of r times an anti-holomorphic function of 
f which, being non singular modular functions, must therefore both reduce to constants 
independent of r, f . The modular integral over the zero modes factor can then be expressed 
by a general formula involving automorphic functions for the moduli fields associated to 
the fixed 2-d torus. For decomposable 6-d tori, one finds [3]: 

A M (T) T ) = ^^ ( ^ )G/ln[m + T i) | r7(Ti) |4 ]) (15) 

G' i=l 

where the sum over G' runs over the distinct N = 2 suborbifolds G' or modular orbits and 
the coefficients denote the associated massless modes beta function slope parameters 

I G f I 

multiplied by the ratios of point groups orders, 1 j g j. The dependence on the Dedekind 
function reflects the target space duality symmetry under the SL(2, Z) modular group. 
The model dependent coefficients can also be represented as [5]: 6 a * = — k a 5 GS , with 
b'a = k[ c (G a ) — J2roc (l+2n l a )c(R a )], where n l a are the massless modes modular weights and 
5 GS the coefficients of the anomaly cancelling Green-Schwarz counterterm. The splitting 
ba = ba + ka5 GS exhibits the characteristic property of the mechanisms responsible for 
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the cancellation of the sigma model duality symmetry anomalies (proportional to 
which involve both threshold corrections (6 a l ) and a gauge group independent Green- 
Schwarz counterterm corresponding to a one loop redefined dilaton field, S + S — > S + 

S + Yli (![t)$ + Ti). For non-decomposable tori, the target space modular symmetry 
is lowered to subgroups of PSL(2, Z). Similar expressions to eq.(15) continue to hold, 
differing by a non-trivial dependence on the sets of allowed moduli, in particular, involving 
rescalings such as % -> T-/3 or T-/4 [28]. 

The moduli independent contributions 5 a are associated with the vanishing of all 
components of the momentum and winding modes, R , yielding therefore a trivial zero 
modes factor equal to unity. No analytic simplification for the modular integral is known 
to exist in this case, for which one must resort to a numerical evaluation. This task is the 
subject of next section and represents the main new result reported in this paper. 
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3. NUMERICAL RESULTS 



Before presenting the results we digress to describe how we deal with the numerical 
integration over the complex parameter r. The two dimensional modular integral can be 
separated in two ways: 



f (frfin,^) = [ 2 dn f 1 dT 2 (f(T U T 2 ) +f(-T 1 ,T 2 )) 

F JO J(1-t?)1 V / 



1 

dT 2 1 dT 1 [f(T 1 ,T 2 )+f(-T 1 ,T 2 )). 
V3/2 JRe(l-T*)? 



(16) 



The general structure of the integrand is that of an infinite sum of terms involving products 
of functions of q, q reading schematically, 

B a (r) = Y J C a (\^x{q)^M = J2 ™a(h L ,h R )q h ^q h * . 

A,/i h L ,h R 

The projection on the modular group invariants is an essential element here in cancelling 
the terms with negative powers of (|) = e ± 27 ™ r i- 27rT 2 j thus leading to non singular expan- 

p2 

sions with powers identified with the conformal weights, hi, = N L + ^ + E — 1, h R = 

Nr + ^ + E q -\, [E = \ EiWCi-M). < \ kv i\ < !]• The functions of r 2 obtained 
upon integration over Ti, as exhibited by the second equation in eq.(16), have discontinu- 
ous derivatives at r 2 = 1, as illustrated in figure 1. When fiL,hn take integral values, the 
r\ (Fourier) integral for t 2 > 1 extends over one period and so selects the terms Kl = h,R. 
The untwisted sector contributions have this property and thus reduce for t 2 > 1 to con- 
stants. The twisted sectors contributions allow (positive) rational values ^ for Iil, h,R and 

2"7T k T2 

so result in decreasing exponentials of the form, e ^~ . Once the constant parts in the 
full integrand, which are identified with the massless modes contributions given by b a , are 
removed, the subtracted integrands (k a B a — b a ) are fastly convergent functions. A cut-off 
at, say, r 2 = 2.2 is more than sufficient to retain the dominant part of the quadrature. 
Nevertheless, the projections involved in the summation over the orbifolds subsectors cause 
strong cancellations which adversely affect the accuracy of final results. The most appro- 
priate way to organize calculations here would be to express analytically the integrand 
in power expansions in q, q prior to the numerical integration [8,9]. However, this proce- 
dure is difficult to implement in a systematic way. We have chosen instead to perform 
all calculations by brute force numerical means and convinced ourselves by various cross 
checks that one could maintain a numerical accuracy better than 10 -2 for orbifolds Zn or 
Zn x Zm, provided that N, M < 6, since the rounding errors grow with the orbifold order. 
The numerical integrations are carried out in the order indicated by the second equation 
in (16). 

Let us quote here useful results concerning the inputs for some of the orbifold param- 
eters. Details regarding the gauge symmetry groups and the massless spectra can be found 
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by consulting refs. [14,24,25]. For the Z 3 j prime orbifolds, the degeneracy factors x(<7,/i) 
count the number of fixed points. Thus, for twisted sectors, x(9i h) = — 27, —9, —3, —1 [g ^ 
1], independently of [h = 1, • • • , 9 N ], for the Z 3 orbifolds with 0, 1, 2, 3 inequivalent Wil- 
son lines, respectively. For the Z7 orbifolds, x{9->h) = ~7(— 1), independently of [h = 
l,---,^], where the first (second) numbers refer to cases without (with) Wilson lines. 
(The reduction of the degeneracy factors in the presence of Wilson lines, reflecting the 
distinguishability of subsets of twisted subsectors, is compensated by a summation over the 
winding numbers, y.) In the Z4 orbifolds, in the absence of Wilson lines, ^t ' 1 ' 2 ' 3 !) = 

-16, x(6> 2 , fll ' 1 ' 2 ' 3 ]) = [16,4, 16,4]. In the Z 3 x Z 3 orbifold with one Wilson line associ- 
ated with the first factor, as in the example presented below, x(g, h) = [3, 3, 3, —9, 3, 3, 3, —9] 
for g = [0 1 ,0l0 2 ,0 1 2 ,0l0 2 ,0l0 1 0l0l0%}, independently of h = 0^0 2 2 . Note that 
x{l,h) = \x(h, 1)| and x(# m , h) = x(@ N ~ m > h.)- The minus signs in the degeneracy fac- 
tors are inserted above in order to account for a twisted sector dependent phase factor 
associated with the chirality. 

The N = 2 subtwisted sectors associated to given (g, h) simultaneously fixed planes, 
consist in the Z 4 orbifold case of a single modular orbit O of (g, h) sectors given by: 
O = {(1, # 2 ), (# 2 , 1), (# 2 , # 2 )}, and in the Z 3 x Z 3 orbifold case of three orbits 0;, associated 
with the three fixed planes, given by: 

1 = {(1, ^ 2 ), (0 2 , 0° 2 > 1 > 2 ), (01 0° 2 > 1 > 2 )}; 2 = O 1 [0 2 - 6 1 ]; 

o 3 = {(i, eAX (i, e\e 2 \ ( (^ 2 ) , i), ( , ^ 2 ), ( , ele 2 )}. 

We present our results for three cases associated with standard embedding (2,2) orb- 
ifolds in Table 1. Results for four non-standard embedding (0, 2) orbifolds are presented in 
Table 2. Details concerning the gauge group and the massless spectra can be found in the 
second reference in [20] and in ref.[5]. Finally, to elucidate the role of discrete Wilson lines, 
threshold corrections results for four realistic cases of orbifolds with three chiral matter 
generations are presented in Table 3. Cases A-C refer to Z 3 orbifolds. Up to extra U(l) 
factors, the observable sector gauge group for Case A [14] coincides with the standard 
model gauge group, while that for Case B, also due to Font et al. [14], is a left-right 
chirally symmetric gauge group extension, SU(3) C x SU(2)l x SU(2)r and that of Case 
C, due to Kim and Kim [29], is an intermediate unification gauge group SU(3) C x SU(3) W . 
Case D in Table 3 refers to a Z 3 x Z 3 [14] orbifold with an observable sector gauge group 
SU(3) C x SU{2) L x SU(2) R x SU{2). 

One of the first calculation of the moduli independent threshold corrections was that 
attempted by Kaplunovsky [1] for the simplest case of standard embedding orbifolds. He 
reported a small gauge group dependent term, A = — A h a Z.^ — 0.07. The Z 3 orbifold case 
with two Wilson lines, designated in Table 3 as Case A, was recently considered by Mayr et 
al., [8]. Assuming tentatively the following decomposition 5 a = —b a A + k a Y, with the first 
component proportional to the factor groups slope parameters and the second to the affine 
levels, they find: A ~ 0.079, Y ~ 4.41. As for the comparison with the existing estimates 
made in fermionic constructions of 4-d superstrings, this is not very teaching because the 



15 



threshold corrections in the models discussed in ref.[6] (A(SU (5)) —A(U (1)) = —24) and in 
ref.[9] (A(SU(3)) — A(U(1)) = —2.5) arise from moduli dependent contributions in N = 2 
sectors only. In a recent systematic study, Dienes and Faraggi [34] report results for several 
new cases. They indicate, in particular, that the above quoted threshold corrections in 
the flipped SU(5) case [9] must be reduced by a factor 3. Let us note here that the models 
obtained in the fermionic construction refer to specific points in the moduli space for which 
one lumps together the moduli dependent and independent contributions. 

The conclusions we draw from Tables 1-3 do not strictly agree with those of refs.[l,6]. 
In our results the component — b a A proportional to the slope parameters is much smaller 
than that quoted above. The coefficient A is never larger than a few % and its sign and 
magnitude change from one group factor to the other and also from case to case. This 
is clearly seen on the corrections 5 a to U(l) factors where the A component is amplified 
by virtue of the larger value taken there by the slope parameters. The analysis of the 
structure of 5 a does not quantitatively support the conjecture made in ref.[8] concerning 
a universal decomposition into two components proportional to b a and k a . We do find, in 
agreement with ref. [8] , a large contribution proportional to the affine levels of approximate 
size, Y ~ 1 — 3. This is not universal, however, but shows rather a tendency to increase 
when including Wilson lines. We remark at this point that the cases in Tables 1 and 2 
featuring significantly enhanced values of S a for certain group factors are precisely those 
cases which involve oscillator states charged with respect to these group factors. Thus, 
charged oscillator states appear as the main responsible for non universal effects. 

We have also examined for the Z 3 x Z s orbifold, the effect of the discrete torsion 
factor, e(TOi,TO2,ni,n 2 ) = e 2 ™^ m ^- m * n ^l N , \p = 0, • • • , N]. The results in Table 3 
refer to the case p = 0. Although the spectrum and hence the slope parameters b a are 
known [14] to depend on the torsion, we find here that the threshold corrections remain 
remarkably stable with variable p > 0. 
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4. UNIFICATION AND ANOMALOUS U(l) SCALE 



4.1 Threshold corrections 



In this section we examine the viability of the perturbative superstring unification 
within the orbifold approach. Let us first discuss the implications of the results obtained in 
Section 3 for the moduli independent threshold corrections. Assuming the simple formula, 
S a = —b a A + k a Y ', then as already noted in connection with eq.(4), one can absorb the 
string threshold corrections into an effective unification scale M' x and an effective string 
coupling constant g' x . Since 5 a are of positive sign, it follows that the moduli independent 
threshold corrections will always result in reduced effective unified coupling constant 
and enhanced (reduced) unification scale, depending on whether the beta function slope 
parameters b a are positive (negative), or equivalently, gauge (matter) dominated. Using 
the numerical values for b a and 5 a in Tables 1-3, we find very small moduli independent 
corrections to the unification scale and/or coupling constant, which attain at most a few 
% . 

Identifying the string moduli independent threshold corrections obtained here, |^ ~ 
0.4, tentatively with a corresponding field theory threshold correction of typical structure 
[30], 8{^) = ±0(1) In H^-, yields for the ratio of the average heavy particle mass to 

unification mass, M H /M X ^ \. Thus, one checks that these contributions are of the same 
order of magnitude as the two loop field theory renormalization corrections [31]. We 
conclude therefore that the moduli independent threshold corrections should mildly affect 
the high energy extrapolation of the gauge coupling constants. More quantitatively, one 
can estimate the corrections to the weak angle and color coupling constant by means of 
the formulas [5], 



sur 6 w (m z ) = - z- + 



ki + k 2 4:7i k\ + k 2 



A , "^7 



-It \ k 
a s (m z ) 



h + k 2 



1 B , ml A B 

— + — log — f- + — - 

a(m z ) 4tt M\ 4tt 



(17) 



where we use the notations: A = — (61^2/^1 — 62), B = — (bi +b 2 — 63 (&i + ^2)/^), A a = 
-{Aik 2 /ki - A 2 ), A B = -(Ai + A 2 - A 3 (k t + k 2 )/k 3 ). Evaluating the threshold 
corrections for Case A in table 3, using k\ = 11/3, yields: 



-4i 



5(sm 2 6 w {m z )) ^ 1.5 x 10" 4 , 8{a- 1 {m z )) ^ 2. x 10" 2 , [5a s (m z ) ~ 2.7 x 10 

where we have set oT x {m z ) = 127.9±0.1. We see that the corrections are rather small and 
lie well within the present experimental uncertainties on these parameters [31], a s (m z ) = 
0.120 ± 0.010, sin 2 6 w {mz) = 0.2324 ± 0.0006. The extreme smallnes of the effect here 
is due to the cancellation of the predominant level dependent component k a Y in S a in the 
linear combinations appearing in A^^. 
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Turning to the moduli dependent corrections Ai m ' ) , we note that these are generi- 
cally of opposite sign with respect to 5 a and so have an opposite effect on the effective 
unification parameters. These contributions become sizeable only to the extent that large 
moduli VEVs and large ratios b' a /b a are used, as is clearly demonstrated on the following 
approximate formula, valid for large VEVs, 



M' x ~ M x 



£ 6 2b n- 



T + T 



(18) 



To estimate the corrections in eqs.(17), one can use the approximate formulas, Aa,b — 

(z',)(\n(2T R ) - |Tr), where (^) = (^-Ib) such that (b) = C% 5 ) for the minimal 
supersymmetric standard model and SA, SB depend on the modular weights parameters 
assignments. The solutions reported in refs. [5,12] give: (^) ~ (2^40)' or equivalently, 
B ,j ~ ( J. In order for these corrections to sin 6w and ot s to reach an order of 
magnitude higher than those found above from the moduli independent corrections, one 
needs at least, Tr = Re(T) = O(10). 



4.2 Standard-like superstring unification scenario 



We shall now present an extended analysis of the string unification picture in which 
the coupling constants scale evolution proceeds through an intermediate threshold at Ma 
induced by an anomalous U(l) mechanism. A two-stage scale evolution is considered: An 
initial short evolution from Ms to Ma, described by the slope parameters set at the 
values predicted in the orbifold models, followed by a wide scale evolution from Ma to mz 
described by the minimal supersymmetric standard model slope parameters. The relevant 
formula reads: 

= M ^ + Y) + 2K In JL + 2bi In £± + A<">(T, f ). (19) 

We regard the five parameters [gx, k\, T , Y = Yq + Y, Ma], which enter explicitly eq.(19), 
as adjustable parameters. Note that Mx has a fixed linear dependence on gx which is 
specified by eq.(3). The compactification scale can be tentatively identified in order of 
magnitude by writing: 

2tt M s fC orb Y 2VC orb M x 



where the compactification radius R and moduli VEV, < T >= T are related as T = 
a°(2^ 2 ' with C or b a calculable constant of order unity [26]. For, say, the Z 3 orbifold, 
C orb = y/S/A. One concludes from eq.(20) that M c /M x ^ l/VT. 

A rough order of magnitude estimate for the anomalous U(l) Higgs mechanism scale 
Ma can be obtained by imposing the condition of a vanishing D-term scalar potential [17], 
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—Da/9 2 a = E a Qll^al 2 + 4a Vfc A ' ^ or a g rou P factor ^7a(1) distinguished by the index A. 
(The triangle anomalies coefficient ca is defined as 48tt 2 ca = Tt(Qa) = 4Tr(Q A ), where 
the traces extend over the massless modes. This enters the Green-Schwarz counterterm 
through the substitution for the dilaton field, S + S — > S + S + caVa, whose function is 
to cancel the various £7^(1) group factor (gauge and gravitational) triangle anomalies, by 
assigning to the gauge vector and dilaton chiral supermultiplet fields the transformation 
laws, Va — > Va — A a — A. A , S — > S + ca^-a-) The predicted magnitude for the scale is: 



if ^ , ^ M P g x 



g x Trace(Q A ) 



192n 2 Q A VkA 



(21) 



Using tentatively for the model dependent ratio the estimate — Tr(Q a) / {QAa\fk~A) — 10, 
one obtains: Ma — l.2g 3 ^ 2 x 10 17 GeV, which indicates that Ma should be of the same 
order of magnitude as M X - 

We use the known experimental values of the gauge coupling constants at the Z-boson 
mass, namely, g\{uiz) = 0.127, g^iynz) = 0.425, g 2 (mz) = 1.44, as inputs to determine 
via eq.(19) three among the above quoted adjustable parameters. We choose these to be: 
gx, Y, Ma- This choice is motivated by the fact that the dependence on these parameters 
in eq.(19) can be made linear by means of an obvious change of variables. The solutions 
for gx,Y , Ma are determined as a function of the remaining free parameters, namely, T 
and fei, and the sets of slope parameters, b^,b' a . For a solution to be acceptable it must 
comply with the perturbation theory constraints that gx and Y be of order unity and with 
the obvious inequalities between scales, < 1, |f^- < 1, which we shall eventually 

supplement by the inequality, < 1, reflecting the assumption that the mechanism 
inducing the scale Ma is a consequence of compactification. 

We shall present the results of numerical applications only for Case A in Table 3, 
setting b a = (—11,-1,3), corresponding to the minimal supersymmetric standard model 
, = (—71.5, —18, —9), as obtained from Table 3, and b' a = (18, 8, 6), 5qs = 7, where the 
choice of slope and Green-Schwarz parameters b' a = & a \ 5gs = X)i ^hs f° r m °duli 
dependent threshold corrections is based on the solutions reported in ref. [5] (see also ref. 
[32]). Regarding k\ as a free parameter when this is predicted to be 11/3 and including 
moduli dependent threshold corrections in a case (such as the Z% orbifold) where these are 
absent, is certainly liable to criticism. However, because the orbifold order appears to have 
a minor influence on threshold corrections and in view of the wide freedom expected in 
the hypercharge gauge coupling constant normalization, we hope that these shortcomings 
do not affect the consistency of our procedure. 

Our main purpose is to explain the non trivial interplay between the various pa- 
rameters which are most significant for string phenomenology. Choosing the particular 
subset, fei, T, as our free parameters while adjusting the others (Y, Ma, gx) to the inputs, 
g a {m 2 z ), [a = 3, 2, 1] is only a technical convenience. Let us first discuss some qualitative 
features of the solutions and, in particular, the correlations among the parameters. The 
dependence on Y and gx shows clearly that any change in gx can be compensated by a 
negative contribution to Y. A decrease of k\ widens the distance between the quantities 
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(fi'a^a) -1 an d so can be compensated by decreasing Ma/M x or g X - Finally, because the 
functional dependence on Ma and gx in eq.(19) involves a logarithm of these quantities, 
one expects a strong sensitivity of the parameters on the inputs. 

The results are displayed in figure 2. These represent a continous two parameters 
(ki,T) family of solutions for gx,Y,MA consistent with a high energy extrapolation 
of the gauge coupling constants joining roughly at the common value, — 25. The 

9 a K a 

physical constraints on Y,gx,MA select a reduced domain for the free parameters, k\ G 
(1.4, 1.8), T E (1,30). The variations with respect to these parameters are monotonic. 
For fixed T, increasing k\ leads to a rapidly (algebraically) increasing Y from large negative 
to positive values and to less rapidly increasing Ma/Mx and gx- Strong variations are also 
found for the T— dependence. However, as T increases past T 25, Y becomes positive 
and nearly independent of k\. The values of k\ on the lower side, k\ < 1.4, are excluded 
by the constraints on Y and those on the higher side, k\ > 1.8, by the constraints on gx 
and Ma/Mx - 

A wide class of solutions occurs with g x « 1 and — Y » 10 3 , independently of T 
and k\. These arise through an obvious compensation effect of the moduli independent 
corrections with gx in eq.(19). Although the Y$ component of Y remains uncalculated so 
far, it appears unlikely that this can much exceed the component Y which was evaluated in 
Section 3 to be of 0(1). In fact, since large Yq is only possible for a strongly coupled string 
theory involving large gx, the above must be regarded as an inconsistent class of solutions. 
(However, because the generic dependence on coupling constant of non perturbative effects 
is expected to be less suppressed in string theory than in field theory [16], e~ c l gx versus 
e -(47r) /g x ^ one cou [ c j poggibly achieve large Y with not too large gx-) In the following 
we shall restrict ourselves to the conventional framework where one assumes a smooth 
connection between string theory and its low energy limit and hence retains the constraints 
gx = 0(l), |Y| = O(10). 

Examining the variation of the solutions with k\ in figures 2(a-c), we see that these 
are very rapid, especially that of Y . The condition Y = 0(1) can be satisfied only through 
a very careful fine-tuning of k± for fixed T, or of T for fixed k\. This is possible only 
in cases where Y changes sign in the relevant intervals of k±,T. The moduli dependent 
corrections are quite essential to achieve a high energy extrapolation consistent with 
superstring unification. Incorporating the threshold Ma provides solutions with reduced 
T. The constraints on Y and gx require 15 < T < 30 and 1.5 < k± < 1.8. Incorporating 
the constraint < 1 restricts this interval to 1.5 < k\ < 1.7. (Narrower intervals would 

be imposed if one also sets lower bounds, say, -j^ > 10 _1 and gx > 10 _1 .) If one takes 

into account the additional constraint < 1, this would lead to the stronger bound, 

~ jj^^/j; < -jf, which would select the narrower interval, 1.5 < ki < 1.6. 

For concreteness, we show in figure 2 (d) the scale evolution of the gauge coupling 
constants for one particular solution as_ determined by the above procedure. One should 
not be disturbed by the large value of \Y\ used here, since the nearby solution determined 
with a carefully tuned value of k\ or T so as to give Y — 0(1), would yield nearly identical 
flows for the gauge coupling constants. This figure illustrates one of the characteristic 
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implications of string unification, namely, that the simultaneous equality at some scale 
of the extrapolated coupling constants has no special significance. The picture depicted 
in figure 2 (d) is rather generic. The most favorable situation corresponds then to an 
approximate joining of the coupling constants flows at a large scale near 5 x 10 16 GeV , which 
is to be identified with the anomalous Ua(1) scale Ma, associated with the decoupling of 
the extra quarks or leptons modes. In the string unification picture, the joining scale 
Ma can be made larger than Mqut because of the slightly reduced normalization of the 
hypercharge group coupling constant and of the spread of the coupling constants at Mx 
which is related to the moduli dependent threshold corrections. 

Let us comment briefly on the sensitivity of the solutions to the slope parameters. 
(Our procedure would obviously break down for b^ ~ b a as this would make the linear 
system of equations, eq. (19), singular) The slope parameters determine the variation 
of the coupling constants from Mx to Ma- The choice of b^ is correlated to that of the 
moduli dependent slope parameters, b' a , since the latter determine the amount by which 
the coupling constants are spread at Mx- Consider first the case of fixed b a . Increasing 
T implies a wider spread of the coupling constants at Mx which should therefore be 
compensated by larger slopes b^ in order to catch up with the extrapolated coupling 
constants up to 10 16 GeV. Rather than showing new plots, we only mention here that if 
one performs a uniform reduction of the slopes by, say, a factor 2, the solutions would 
rule out the entire domain in k\,T except for a narrow region around T = 15, k\ = 1.7. 
Conversely, enhancing the slopes b^ by, say, a factor 2 ameliorates the initial picture 
without changing qualitatively the character of solutions. One concludes therefore that 
the cases involving negative slope parameters b^ with large absolute values (richer matter 
spectra), which are generic in orbifolds model building, are more favorable for unification. 

The choice of b a = b' a — k a dcs is also quite sensitive. Rather than performing an 
exhaustive study we have considered two other cases obtained from ref. [5] and further 
motivated in ref. [32]. Applying the above procedure of solution for these cases, we found a 
significantly worsened picture. The first case, characterized by b' a = (7.5, 2.5, 1.50), Sqs = 
2.5, admits solutions only for large values of T > 30 and correspondingly large k\ > 
1.8. It improves slightly if reduced values are used for the slopes b^- The second case, 
characterized by b' a = (—4.67,4, 5), 5qs = 6, admits no solutions at all, mainly on account 
of an incompatibity between the constraints on Y and Ma/M X - One concludes therefore 
that negative or small values for the N = 2 slope parameters b' a do not constitute a 
favorable option. 

Having focused so far on standard-like compactification models, we briefly discuss 
the other two possible classes of superstring models. The first refers to compactification 
models with grand unified groups, SU(5) [6] or 50(10) [21] (up to extra U(l) factors), 
with a flipped assignment for the matter fields with respect to the standard GUT basis or 
with a regular GUT assignment involving higher affine levels, k > 1 [33]. A perturbative 
weak coupling scenario assuming a smooth evolution from Mqut to Mx can be carried 
out in the manner described above either by setting the parameters, bo, b' G and Y at values 
specified by the models or by imposing appropriate constraints on them. It should not be 
difficult to obtain satisfactory solutions for gx and Ma by following a procedure similar 
to that used above. An alternative strong coupling scenario could also be envisaged [21] 
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if the slope be takes a large (gauge dominated) positive value and gx is large so as to 
lead to GUT group G with renormalization group invariant scale comparable to the string 
scale, A G = M x e~ 87T k G/ b Gg' x . Although such a scenario forbids a smooth connection from 
string theory to the low energy field theory, it still provides a prediction for the GUT scale, 
namely, Mqut — A G . 

The second class of compactification models corresponds to intermediate unification 
on a semisimple electroweak gauge group. One interesting example is Case C in Table 3 
where the gauge symmetry at compactification, SU(3) C x SU(3) W x U(l)p 3 , breaks down 
to the standard model group at an anomalous U(l) scale according to SU(3) W x U(l) p a — > 
SU(2) W x U(1)y, where Y = Ts w + Using the information supplied in ref.[29], we find 
a level parameter k(Ps) = |. This implies a normalization of the hypercharge coupling 
constant such that k\ = 1 + ^ = ||. Although this falls well below the favorable interval 
of k\ values specified above, it is nevertheless interesting that the situation for Case C is 
exactly opposite to that found above for Case A. A simple argument was recently made 
[34] that for any (orbifold or fermionic) model realizing a direct compactification to the 
standard model group, requiring a correctly normalized hypercharge imposes the bound: 
k\ > |. Therefore, an intermediate unification of hypercharge in a non-abelian group 
factor would appear as one of the most viable options. 



5. CONCLUSIONS 



Our results indicate that the moduli independent threshold corrections are compara- 
ble in size to those for gauge field theories in spite of the fact that infinitely many massive 
states are integrated out for superstrings. The corrections are marginally relevant at the 
current precision levels for the low energy gauge coupling constants. The largest contri- 
butions reside in a group independent component k a Y of size Y ~ 1 ~ 3 which remains 
relatively stable with respect to the orbifold order or to the choice of gauge group em- 
bedding and Wilson lines. The component — b a A is much smaller, |A| < 10 -2 and model 
dependent. 

In order for the large value of the predicted string unification scale Mx not to conflict 
with observations, one needs both moduli dependent threshold corrections (with associ- 
ated compactification scale -j^ ~ -7= « 0.3) as well as a weak hypercharge group level 
parameter varying in the narrow interval, k\ = 1.4 ~ 1.7. The information that the mod- 
uli independent corrections are 0(1) is useful in providing stronger correlations among the 
parameters relevant to string phenomenology. Postulating an anomalous U(l) mechanism 
at a scale 0.1 < Ma/M x < 1 significantly eases the above constraints on slope parameters 
while raising the bound on the allowed values of Mq- The resulting picture is intermediate 
between a delayed joining of the coupling constants flows, due to the smaller value of ki, 
and of a continued flow beyond crossing, consistent with the moduli dependent threshold 
corrections. Our analysis emphasizes the need of constructing orbifold models combin- 
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ing the property of a low value for the hypercharge group level parameter along with the 
usual desirable features, namely, three chiral families, low rank gauge group and N = 2 
subsectors. 
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TABLES CAPTIONS 



Table 1. Threshold corrections for the .£3,4,7 orbifolds with standard gauge embeddings. 
The entries in the first line are the rotation angles 9{, [i = 1,2,3] and the shift vectors 
v\[i = 1,2,3] V 1 , [I = l,---,8]. The second and subsequent columns correspond to 
the gauge group factors in the observable and hidden (primed) sectors. For each column 
beyond the first, the first line entry gives the levels k a , the second line gives the beta 
function slope parameters b a or, for the non-prime orbifolds with N = 2 suborbifolds, the 
pairs (6^ =1 ,6' a ), such that b a = b^ =1 + b' a . The third line gives the moduli independent 
threshold corrections 5 a . 

Table 2. Threshold corrections for orbifolds Z3 4 with non standard gauge embeddings. 
For each case, the first line gives the shift vectors V 1 , V 1 , [I = 1, • • • , 8]. The second and 
subsequent columns correspond to a selection of the gauge group factors in the observable 
and hidden (primed) sectors. For each column beyond the first, the first line entry gives the 
levels k a , the second line gives the beta function slope parameters b a or, for the non- prime 
orbifolds with N = 2 suborbifolds, the pairs (6^ =1 ,6' a ), such that b a = b^ =1 + b' a . The 
third line gives the moduli independent threshold corrections 5 a . 

Table 3. Threshold corrections for a selection of three-generations orbifold models with 
two Wilson lines (Cases A-C) and one Wilson line (Case D). For Z 3 orbifolds, the winding 
number parameters attached to the Wilson lines take the values: rrikj = 0, ±1. Case A 
is a standard model group Z3 orbifold model from Font et al., [14] (section 4.2): 3V 1 = 
(l 4 2000)(20 7 )',3a( 2 = (0 7 2)(0110 5 )', 3a 3 , 4 = (11121011)(110 6 )'. Case B is a left-right 
group Z 3 orbifold model from Font et al., [14] (section 4.3): 3V 1 = (1 4 2000)(20 7 )', 3af = 
(0 7 2)(00110 4 )', 3c4 = (1 3 21 3 0)(110 6 )'. Case C is an intermediate unification group Z 3 orb- 
ifold model from Kim and Kim [29]: 3V 1 = (11211200) (0 8 )', 3a{ = (0 3 11211)(1 4 4 )', 3a^ = 
(0 7 2)(1 8 )'. Case D is an intermediate unification group Z3 x Z3 orbifold model with 
one Wilson line from Font et al., [14] (section 5): 3v\ = (1, 0, -1), 3*4(0, 1, -1); 3V"/ = 
(2110 5 )(110 6 )',3V/ = (020 6 )(0 - 1111000)'; 3aS 1)J = (0 5 11 - 2)(0 5 11 - 2)'. (The indices 
1, 2 refer to the two Z^ factors.) For each column, the first line entry gives the levels k a , 
the second line gives the beta function slope parameters b a or, as in Case D, the pairs 

( °r, ) such that b a = 6^ =1 + b' a . The third line gives the moduli independent threshold 

a 

corrections 5 a . 



FIGURES CAPTIONS 



Figure 1. The threshold function — S a (r), integrated over n, is plotted as a function 
of T2 for the Z3 orbifold model group factor SU(3) C of Case B in Table 3. We show the 
contributions of the untwisted (continuous line) and of the twisted sectors (double-dashes). 
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Figure 2. One loop renormalization group analysis of superstring unification parameters 
based on high energy extrapolation of the gauge coupling constants starting from their 
experimental values at mz- The solutions for —Y (figure a), Ma/Mx (figure b) and gx 
(figure c) are plotted as a function of k± for a discrete set of values of the moduli VEV, 
T = 1 (continuous), 10 (long dash short double-dashes), 15 (long dash short dash), 20 
(dash dot), 30 (dash). The slopes discontinuities exibited by Y in figure (a) arise because 
of the changes of sign of Y in this semilogarithmic plot. (For the T = 30 curve, Y > 0.) 
We display in figure (d) graphs of the gauge coupling constants (-ff- , [a = 3,2,1]) 
variation with renormalization scale for the particular solution characterized by the values, 
k\ = 1.6, T = 20, yielding the solution Y = —114, M A /M X = 0.38, g x = 0.63. 



27 



TABLE 1 
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TABLE 2 
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TABLE 3 
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